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It is shown the analysis [l] for QED in 2+1 dimensions with N four-component 
fermions in the leading and next-to-leading orders of the 1/N expansion. As it 
was demonstrated in [1] the range of the admissible values N, where the dynamical 

o . 

CN , fermion mass exist, decreases strongly with the increasing of the gauge charge. So, in 

Qh' Landau gauge the dynamical chiral symmetry breaking appears for N < 3.78, that 



< 

On 
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is very close to the results of the leading order and in Feynman gauge dynamical 
mass is completely absent. 



Quantum Electrodynamics in 2+1 dimensions (QED3) has acquired increasing attention 
Q] because of its similarities to (3+1) dimensional QCD. Moreover, last years a new strong 



interest comes to QED3 in the relation with graphene properties (see [8j and discussions and 



references therein). Graphene, a one-atom-thick layer of graphite, is a remarkable system 

00 ■ n 

00 ' with many unusual properties that was fabricated for the first five years ago |9| . Theoretically 



it was shown long time ago [10[ that quasiparticle excitations in graphene are described by 
the massless Dirac equation in (2+1) dimension. This explains why the bilayer graphene in 
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external fields is a subject of intensive recent study 

A number of investigations have been performed for the study of dynamical chiral sym- 
metry breaking in QED 3 and very different results have been obtained. Using the leading 
order (LO) in the 1/N expansion of the Schwinger-Dyson (SD) equation, Appelquist et al. 
(2! showed that the theory exhibits a critical behavior as the number A" of fermion flavors 
approaches A" c = 32/ V 2 ; that is, a fermion mass is dynamically generated only for A" < N c . 
On the contrary, Pennington and collaborators [3j, adopting a more general non-perturbative 
approach to the SD equations, found that the dynamically generated fermion mass decreases 
exponentially with N, vanishing only as A" — > 00. This conclusion was supported also by Pis- 
arski [4] by the use of the other methods. On the other hand, an alternative non-perturbative 
study by Atkinson et al. 1^] suggested that chiral symmetry is unbroken at sufficiently large 
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N. The theory has also been simulated on the lattice j^, 7\. Remarkably, the conclusions of 
Ref. are in the agreement with the existence of a critical iV as predicted in the analysis 
of Ref. Q while the second paper [zj] contains the opposite results. 

Because the critical value N c is not large, the contribution of the higher orders in the 
1 /N expansion can be essential and may lead to better understanding of the problem. The 
purpose of this work is to consider the 1/N correction Q to LO result 

1. The Lagrangian of massless QED 3 with N flavors is 

1 



L = V(id - eA)^ - -F. 
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where ^ is taken to be a four component complex spinor. In massless case, which we are 
considering, the model contains infrared divergences, which can be canceled when the model 
is analyzed in a 1/N expansion 13j. Since the theory is massless, the mass scale is the 
dimensional coupling constant a = Ne 2 /8 which is kept fixed as N — >• oo. 

n 

Following [2J we study the solution of the SD equation. The inverse fermion propagator 
has the form 

S-\p) = -[1 + A{p)\ + 

where A(p) is the wave- function renormalizable coefficient and is a dynamical, parity- 
conserving mass taken to be the same for all the fermions. 
The SD equation is 

2a r d 3 k rD,u(p-k)[l + A(k)]Ck + nk))T^p,k) 

N lr J (2tt)3 [l + A(k)] 2 (ki + ^(k)) ' { > 

where 1 

^ [p) p 2 [i + n(p)] 

is the photon propagator and T u (p, k) is the vertex function. 



2. The LO approximations in the 1/N expansion are 

A(p) = 0, Tl(p) = a/\p\, T»(p,k) = Y 



Following 12[ we introduce a nonlocal gauge-fixing term. The detailed analysis of this possibility has 
been given in Ref. 14 1 . 
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where we neglect the fermion mass in the calculation of II (p). The gap equation is 

E(V) = 8a(2 + Tr f ^ ^ (2) 

KP) N J (2tt) 3 k 2 [(p - k) 2 + a | p - k ]]' K > 

where we ignore the term T?(k) in the denominator of r.h.s. . 
Following [2|, we set 

m = (k 2 T- (3) 

One can see, that for large a the r.h.s. of (T5]) together with condition ([3]) (and the contribu- 
tions of higher orders) can be calculated by the standard rules for massless diagrams of the 
perturbation theory (see, for example, [lsj] ) . Thus, we have for large a 

with (3 = (-a) (a + 1/2) and L = n 2 N, or 

a ± = (-l±[l-16(2 + 0/£] 1/2 )/4. (5) 



We reproduce the solution given n Ref. [2j. That analysis yields a critical number of 
fermions N c = 16(2+0 /tt 2 w 1.62(2+0 (i.e. L c = 16(2+0), such that for N > N c S(p) = 
and 

S(0) ~ exp[-27r/(N/N c - 1) 1/2 ] 



for < N c . Thus, chiral-symmetry breaks when a becomes complex, that is 
3. The next-to-leading order (NLO) approximation has been included in [l|, 



or N < N c . 
using; the 



differential equation method [121]. The results have a cumbersome form [161 ] . which is similar 
to results for complicated massless diagrams obtained using Gegenbauer polynomials 18]. 
In l| we have analyzed simplified form, which contains only the terms ~ (— a)~ k and ~ 
(a + l/2)~ k (k = 1, 2, 3) from the series given in [if]]. These terms are most important in 
the neighborhood of the critical point N c . The eq. (j3J) is replaced now by 

where /«) = 4(1 - - <p(i) = 176/9 - 4tt 2 - (16/3)? + 4? 2 . 

Let us get the exact critical value N c from eq. Supposing a = a c = — 1/4 we obtain 
the critical values in the following form 

N«± = ^[(2 + ± ((2 + 2 + 4/(0 + V(0/4) V2 ], (7) 
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i.e. 



JV Ci+ (0 = (3.31, 3.35, 3.09, 2.81), iV c _(£) = (-0.07, 0.38, 1.29, 1.88) 

for f = (0.0, 0.3, 0.7, 0.9), respectively. 

Notice the intriguing fact that follows from (JTj). The addition of 1/N correction leads to 
the occurrence of the second critical point (for 0.05 < £ < 0.95) such that for N < A Ci _ 
the chiral symmetry does not break. The dynamical mass generation exists in the interval 
between the critical points N c _ and iV c +. For £ > 0.95 this interval disappears and the 
chiral symmetry breaking is absent. For small values of gauge parameter £ (£ < 0.05) new 
critical point does not occur. 

The solution of the eq.flH]) is 

& = ± [2 + £ + ^ ± ((2 + 2 + 4/(£) + 2(2 + 0^ + ^) 1A 



has the simple form in Landau gauge 

1 r. <p(0) 



P±(Z = o) 



L 



l + ^±VT/3 1 + 



3 <p(oy 

14 L 



1 + 



3 y(Q) 
14 L 



(8) 



where the last term in r.h.s. of eq.(jH]) is very small for L ~ L c . Leaving it out we get the 
following equation for (3 + 

Me = o) - 1 + + (i + vwo)) ^ - 2 -f (i - ^) , 



which has coefficients are close to those from the paper [12] . 



Resume. We reviewed the results of [l|, where the 0(1/A^ 2 ) terms have been in- 
cluded into SD equation and the strong gauge dependence of the result has been found. 
Hence, the addition of 1/N correction does not lead to the essential improvement in the 
understanding of dynamical chiral symmetry breaking. However, as it was shown in [l|, 
in the Landau gauge the inclusion of 0(1/N 2 ) terms slightly changes only quantitative 
(but not qualitative) properties of the LO results. Thus, in the Landau gauge the analy- 
sis in [1] gives further evidence in favor of the solution has been given by Appelquist et al. (2]. 
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